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Landau gauge ghost and gluon propagators in SU(2) lattice gauge theory: 
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We reinvestigate the problem of Gribov ambiguities within the Landau (or Lorentz) gauge for 
the ghost and gluon propagators in pure SU(2) lattice gauge theory. We make use of the full 
symmetry group of the action taking into account large, i.e. non-periodic Z(2) gauge transformations 
leaving lattice plaquettes invariant. Enlarging in this way the gauge orbits for any given gauge field 
configuration the Landau gauge can be fixed at higher local extrema of the gauge functional in 
comparison with standard (overrelaxation) techniques. This has a clearly visible effect not only 
for the ghost propagator at small momenta but also for the gluon propagator, in contrast to the 
common belief. 
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I. INTRODUCTION 

It is an important task to compute (Landau 
or Coulomb) gauge gluon, ghost, fermion propa- 
gators and the basic vertex functions from non- 
perturbative approaches to SU(N) gauge theo- 
ries, like Dyson- Schwinger equations or the lat- 
tice formulation. On one hand one is interested 
in their behavior in the infrared limit in order 
to extract non-perturbativc informations on var- 
ious observables, e.g. the QCD running couplings 
a s (q 2 ), to understand quark and gluon con- 
finement within the Gribov-Zwanziger scenario 
d E3- Q3| - or to check the Kugo-Ojima confinement 
criterion for the absence of colored states |j] . On 
the other hand it is technically important to see 
to what extent these different non-perturbative 
approaches provide results consistent with each 
other in the non-perturbative region, i.e. at low 
momenta. At present we are still far from draw- 
ing final conclusions in this respect. In particular 
the Dyson-Schwinger approach 0, j always 
relying on a truncated set of equations, provides 
results which look quite different in the infinite 
volume limit compared with those obtained on 
a torus |3, |l0| , while the latter show at least 
qualitative agreement with recent results of nu- 
merical lattice simulations [Til ]. 

It is well known that gauge fixing in the non- 



perturbative range is faced with the Gribov am- 
biguity problem, which means that there can be 
many gauge copies for a given gauge field satis- 
fying the Landau gauge condition d^A^ = 
within the Gribov region, the latter defined by 
the positivity of the Landau gauge Faddeev- 
Popov operator. In recent years one has checked 
in greater detail how strong Gribov copies can 
influence the infrared behavior especially of the 
gluon and ghost propagators. Several groups of 
authors came to the conclusion that while there 
is a clearly visible influence on the ghost prop- 
agator, the gluon p ropagator seems only weakly 
affected [Hill M, M, El 

Recently Zwanziger has argued that in the infi- 
nite volume limit the influence of Gribov copies: 
"... might be negligible, i.e. all averages taken 
over the Gribov region should become equal 
to averages over the fundamental modular re- 
gion" 3]. However, in practical lattice simula- 
tions we are always restricted to finite volumes. 
Thus, Gribov copies have to be taken into ac- 
count properly before extrapolating to the in- 
frared and infinite volume. 

In this paper we present a reinvestigation of 
the Gribov copy problem for the SU(2) case. The 
usual way to fix the (Landau) gauge on the lattice 
is to simulate the path integral in its gauge in- 
variant form. Subsequently each of the produced 
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lattice gauge fields U = {U x ^} is subjected 
to an iterative procedure maximizing the gauge 
functional 

U°, p = 9(x) U^ g^x + fi) (1) 

with respect to local gauge transformations g = 
{g(x) e SU(2)}. V — L d denotes the number 
of lattice sites in d = 4 dimensions. The local 
maxima of F(g) satisfy the differential lattice 
Landau gauge transversality condition 

(M»)(s) = A%{x + A/2) - A^x - A/2) = , 

(2) 

where the lattice gauge potentials are 

A fi (x + ji/2) = ±- i (U X}ll -Ulj. (3) 

The standard procedure assumes periodic gauge 
transformations and employs the overrelaxation 
algorithm. In what follows we shall abbreviate it 
by SOR. The influence of Gribov copies can be 
easily studied by taking various initial random 
gauge copies of the gauge field configurations be- 
fore subjecting them to the SOR algorithm. 

At this place it is worth to note that the widely 
- at least till now - accepted approach to compute 
e.g. a gauge-variant propagator G is to choose 
always the gauge copy with the highest value of 
local maxima F max (or the best copy) found for 
the gauge functional Q). One can then hope to 
have found a copy belonging to the so-called fun- 
damental modular region or at least being not 
far from it. In order to find the best copy for 
each thermalized gauge field configuration one 
needs to compare ,F maa ;-values for a pretty large 
amount of gauge copies, which is a rather time 
consuming procedure. A reasonable question is 
if the use of only one gauge copy (the first copy) 
provides us with the same - within errorbars - 
values of the propagator as the use of the best 
copy. This logic brings us to compare the propa- 
gator calculated on best copies (G( bc )) with that 
on the first copies (G^ c ^). The relative devia- 
tion SG = |(G (/C ) - G( bc) )/G (bc )| provides then 
a useful quantitive measure of the Gribov am- 
biguity of the quantity under consideration. We 
shall discuss this measure throughout the present 
paper. 

Of course, one can enhance the effect of Gri- 
bov copies by comparing instead the best copies 
with the worst copies, i.e. with those having the 



smallest values F max found from the repeated use 
of a given maximization method. This attitude 
has been taken in Ref. 0] in order to highlight 
a Gribov copy effect for the gluon propagator. 

In Refs. for SU{2) and 11] for SU{3) some 
of us already have thoroughly discussed the im- 
pact of Gribov copies within the SOR framework 
by comparing first and best copies. From this 
point of view the gluon propagator did not de- 
pend on the copies within the statistical noise, 
whereas the ghost propagator clearly was de- 
pending on them in the infrared. But the data 
for the ghost propagator obtained for different 
lattice sizes showed an indication for a weakening 
of the dependence on the choice of Gribov copies 
for increasing lattice size at fixed momentum, in 
agreement with Zwanziger's claim 

Here we enlarge the class of possible gauge 
transformations by taking into account also non- 
periodic center gauge transformations. This will 
allow us to maximize further the gauge functional 
and to see a quite strong Gribov copy effect also 
for the gluon propagator at finite (lattice) vol- 
umes. 

In Sec.^we shall explain the improved gauge 
fixing procedure. In Sec. 11111 we define the prop- 
agators to be calculated. In Sec. lIVI we are going 
to present our results for the gluon and ghost 
propagators, whereas in Sec. the conclusions 
will be drawn. 



II. IMPROVED GAUGE FIXING 

We shall deal all the time with SU{2) pure 
gauge lattice fields in four Euclidean dimensions 
produced by means of Monte Carlo simulations 
with the standard Wilson plaquette action. We 
restrict ourselves to the confinement phase at 
T = 0. 

To fix the gauge we employ the standard Los 
Alamos type overrelaxation with uj = 1.7. 

Our generalization of the standard gauge fix- 
ing procedure SOR comes from the simple obser- 
vation that gauge covariance for periodic SU(2) 
gauge fields on a d-dimensional torus of exten- 
sion L d allows gauge transformations which are 
not necessarily periodic but can differ by a group 
center element at the boundary: 

g(x + L£>) = z l/ g(x),z v = ±leZ(2). (4) 

In light of this it is legitimate to allow, during 
the maximization of the gauge functional in the 
gauge fixing procedure, for gauge tranformations 
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which differ by a sign when winding around a 
boundary. Let v be the direction of such bound- 
ary. Any such gauge transformation can be de- 
composed into a standard periodic gauge trans- 
formation (which we may call a "small" one) and 
a flip of all links U v (x) — > — U v (x) of a 3- 
plane at a given fixed x v . Given a "small" ran- 
dom gauge copy of the configuration we have thus 
performed a pre-conditioning step for the gauge 
functional by sweeping in every direction all 3- 
planes in succession and comparing the value of 
the flipped with the unflipped gauge functional. 
The flip is accepted if the gauge functional in- 
creases. It is easy to see that such a procedure 
is independent of the order of choosing the 3- 
planes and that only one sweep through the lat- 
tice is required to maximize the functional. The 
gauge copy obtained at the end of this procedure 
is then used as a starting point for a standard 
maximization procedure. We call the whole pro- 
cedure FOR. 

Analogously to the SOR method the FOR pro- 
cedure can be repeated with different initial ran- 
dom gauges in order to find a best copy ( be ) 
in comparison e.g. with the first random copy 
( f c ) . We shall check the convergence of the 
be -propagator results for the best copies as a 
function of the number n copy of random initial 
copies. 

III. GLUON AND GHOST 
PROPAGATORS 

We turn now to the computation of the gauge 
variant gluon and ghost propagators within the 
Landau gauge. 

The lattice gluon propagator D^Kp) is taken 
as the Fourier transform of the gluon two-point 
function, i.e. the expectation value 

DfM = (A;(k)A h l/ (-k)) u (5) 

A a ^(k) is the Fourier transform of the lattice 
gauge potential A^(x + fi/2). p denotes the 
four-momentum 




with the integer-valued lattice momentum k^ S 
(-L/2, +L/2]. a is the lattice spacing. 



The lattice ghost propagator is defined by in- 
verting the Faddeev-Popov (F-P) operator, the 
latter being the Hessian of the gauge functional 
Eq. (|TJ. The F-P operator can be written in 
terms of the (gauge-fixed) link variables U x<ll 
as 

n rah \ , Aab s r>ab r /~mb x 

xy ~ °xy~ °x+fi,y — '-'x,ii °x~)J,,y 

(7) 

with 

A-x b ,n — 2 ^ ab ^ [^ x ^ ^ic— A,/t] ' 

B£» = \ Tr [aV U x>fl ] , 

Cll = \^[a a a b U x ^] , 

where the a a , a = 1,2,3 are the Pauli matrices. 
In the continuum M c x lb corresponds to the op- 
erator M ab = -d^Df, with D ab the covariant 
derivative in the adjoint representation. 

The ghost propagator in momentum space is 
calculated from the ensemble average 

x,y 

= 8 ab G( P ) . (9) 

Following Ref. 0, we have used the conju- 
gate gradient (CG) algorithm to invert M on a 
plane wave ip c = { S ac exp(2ni k-x) }. 

After solving M<fi = ip c the resulting vector 
is projected back on ip so that the average 
G cc (p) over the color index c can be taken ex- 
plicitly. Since the F-P operator M is zero if 
acting on constant modes, only k ^ (0,0,0,0) 
is permitted. Due to high computational require- 
ments to invert the F-P operator for each k, sep- 
arately, the estimators on a single, gauge-fixed 
configuration are evaluated only for a preselected 
set of momenta k. 



IV. RESULTS 

We consider various bare couplings in the in- 
terval (3 = 4/<7o €[2.1,2.5] and lattice sizes up 
to 20 4 . We compare the gluon and ghost propa- 
gators obtained with the alternative gauge fixing 
methods SOR ('flips off) and FOR ('flips on') 
both for the first (fc) and best copy (be). In 
order to find the best copies we always generate 
20 initial random gauge copies. 
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FIG. 1: Gluon propagator (left) and ghost propagator (right) at lowest momentum p m in ~ (2 / a) sm(n / L) 
versus number of random copies employing the FOR method ('flips on') at /3 = 2.5 and 2.4, respectively 
(lattice size 16 4 ). 



In Fig.^ wc illustrate for the FOR method how 
fast the gluon and ghost propagators are converg- 
ing when determined from the best copy out of 
the first n copy copies. We see plateaus occur- 



ing for 



> O(10). We have convinced our- 



selves that O(20) copies are sufficient at least for 
/3 > 2.3 and lattice sizes up to 20 4 . For the SOR 
method the convergence is faster - although to 
worse values of the gauge functional - such that 
in principle a smaller number of copies would be 
sufficient within the given parameter range. 

Mostly we have concentrated on the low- 
est non-trivial on-axis lattice momentum 
Pmin — (2 / a) sm(n / L) and some multiple 
on-axis momenta in order to study the infrared 
limit for given lattice size and bare coupling. We 
are aware of the fact that this choice is by far too 
restrictive in order to get reliable results for the 
(renormalized) propagators in the continuum 
and thermodynamic limit. 

In Fig. |21 and Fig. [3] we show our results for 
the lattice gluon D(p m i n ) and ghost propagators 
G(pmm) for 12 4 and 16 4 lattices, always for the 
smallest non-vanishing momentum. In order to 
demonstrate the effect of the Z(2) flips in com- 
parison with the SOR results obtained with be 
and fc copies we show three sets of data 
points: black dots correspond to FOR - 'flips on' 
and be copies and open circles (squares) corre- 
spond to SOR - 'flips off' for be (fc) copies. 



The corresponding data are listed in Table |U 

We clearly see that the FOR method leads 
to an additional visible Gribov copy effect not 
only for the ghost propagator but also for the 
gluon propagator. The effect is even more pro- 
nounced at higher (3- values, i.e. at smaller 'phys- 
ical' lattice sizes. We have convinced ourselves 
that this is compatible with the behavior of the 
average maximal gauge functional (F max ). Its 
relative difference determined with be copies for 
the FOR method versus the SOR method is also 
rising with j3. Later on we shall see that this 
observation is also in one-to-one correspondence 
with the gauge copy dependence for fixed j3 and 
varying lattice size. The anatomy of the (new) 
FOR gauge copies deserves further studies in the 
future. 

In order to illustrate the strong Gribov copy 
effect in a slightly different manner we compare 
smoothed distributions for the mean value esti- 
mators for the gluon and ghost propagators for 
the be with the FOR and SOR method, re- 
spectively (see Fig. |3J. The mean value distri- 
butions have been obtained in accordance with 
the bootstrap method from replica of se- 
quences of randomly selected data. Such boot- 
strapped resampling was applied to the initial 
MC data set as a whole, the amount of repli- 
cas being typically 200. To smoothen the dis- 
tribution we have used the standard Nadaraya- 
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FIG. 2: Gluon propagator D(pmin) at lowest momentum for various j3 and for lattice sizes 12 4 (left) and 
16 4 (right). Full dots refer to FOR f c and open squares (circles) correspond to SOR fc (be). 
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FIG. 3: Ghost propagator G(pmin) as for Fig. [5] 



Watson method with normal kernel and an 
improved Silverman's rule of thumb for the choice 
of the corresponding bandwidth. 

It is worth mentioning that the statistical er- 
rors for most of our data have also been estimated 
through bootstrapped resampling. 

We have also studied how the Gribov copy ef- 
fect develops for larger momenta p(k). We have 
used multiples of the minimal lattice momentum 



k = (0, 0, 0, fc), k — 1, 2, 3, 4 along one axis. We 
compare for the gluon propagator the be SOR 
results with be FOR results in terms of the rel- 
ative deviation 

SD( P ) = (D^ R - I^Mok - (10) 

and analogously for the ghost propagator G(p) 
at various /3-values and with fixed lattice size 
16 4 (see Fig- EJ- F° r the gluon propagator our 
results are restricted to only one (3- value because 
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TABLE I: Data for the gluon propagator D(p) (left) as well as for the ghost propagator G(p) (right) at lowest 
momentum p — pmin obtained with FOR (be) and SOR (be and fc) methods on 12 4 and 16 4 lattices. 




26 27 28 29 30 31 32 



FIG. 4: SOR- and FOR-distributions for D {bc) (p min ) 

of the much stronger statistical noise. Neverthe- 
less, the results presented for the gluon propaga- 
tor point into the same direction as for the ghost 
propagator. The effect of Gribov copies still re- 
mains noticable at p > p m in, although decreasing 
for rising momenta. The data for the ghost prop- 
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eft) and G (bc) (p m i n ) (right) at P = 2.5 and 16 4 lattice. 

agator at various momenta obtained from inde- 
pendent Monte Carlo runs are also collected in 
Table HIl 

We have also made a corresponding check for 
the gluon propagator at zero momentum. On a 
lattice of size 20 4 and for the same (3 = 2.5 we 
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FIG. 5: Left: relative deviation SD(p) = (-DgoR — ^toh)/-^toh m percent for the gluon propagator at 
various (on-axis) lattice momenta p(k) (lattice size 16 4 , P = 2.5). 

Right: the analogous relative deviation for the ghost propagator for the same lattice size but for /3 = 2.2, 2.3, 2.4 
and 2.5. 



observed a deviation between the be FOR and 
SOR results of the order 0(25%). This would 
of course have consequences for estimates like in 
Ref. [l^ . |20| , since the infinite volume extrapola- 
tion of -D(O) there performed, although probably 
remaining finite, will definitely suffer from un- 
controlled systematic uncertainties. 

It is interesting to study the volume depen- 
dence of the Gribov copy effect, in view of 
Zwanziger's recent claim mentioned at the be- 
ginning Q . First of all we have convinced our- 
selves that the number of gauge copies is strongly 
rising with the lattice volume as it should be. 
This is clearly demonstrated in Fig. providing 
the distributions of the number of gauge copies 
per configuration found with the FOR method 
('flips on') for lattice sizes 8 4 and 16 4 at j3 = 2.40. 
In both cases we have generated 100 configura- 
tions with 100 gauge copies each. It turns out 
that identical (or degenerated) copies can be well 
recognized at an accuracy for the gauge func- 
tional Eq. (JTJ of O(10 _1 °). Adjacent copies nor- 
mally differ in the values for the gauge functional 
at a level of O(10 -6 ). Now let us compare the 
distributions of the corresponding values of the 
functional F for each copy found. In order to 
normalize the values with respect to the highest 
(i.e. best) value per configuration we show the 

relative deviation (F^al -F max )/F, 



5 10 



b 5 





8 4 lattice; p=2.4 




Flips ON 




n N„ nr =100; N cop =100 






E=10-'° 


rf 






in 



R 10 



b 5 




4 b a l Thefre- 



10 20 30 40 50 60 70 80 90 100 
16* lattice; (5=2.4 

Flips ON - 
N conl =100; N com =100 
E =10-'° 

r nn.n.nrm 

10 20 30 40 50 60 70 80 90 100 
# of different copies per configuration 

FIG. 6: Distributions of the number of different 
gauge copies found with the FOR method at /3 = 2.40 
for lattice sizes 8 4 and 16 4 . 



quency distributions of these values are shown in 
Fig. 13 for the same ensembles as used for Fig. 
There is a very clear tendency that the variance 
of the gauge functional becomes much smaller 
if we increase the lattice volume. A similar ten- 
dency becomes visible in Fig. [HI where we plot for 
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15.6(1) 


2.87(1) 


1.171(1) 


0.6690(3) 






be 


14.8(1) 


2.83(1) 


1.165(1) 


0.6673(3) 


2.50 


200 


fc 


14.1(2) 


2.586(8) 


1.090(1) 


0.6359(4) 






be 


13.4(1) 


2.564(6) 


1.088(1) 


0.6352(3) 



TABLE II: Ghost propagators G(p) on the 16 4 lattice 
for various on-axis lattice momenta p(k). 



the same set of configurations and gauge copies 
the distributions for the single values of the ghost 
propagator for the lowest non-vanishing on-axis 
momentum. Also in this case we have normal- 
ized the single values as (G^ - G)/G {bc \ i.e. 
taking the relative deviation of the propagator 
at a given copy G from the value computed on 
the best copy G bc , the latter chosen again with 
respect to the gauge functional value. We see 
that the long tail seen for the smaller lattice dis- 
appears for the larger lattice. Although the fact 
that close values of the gauge functional will not 
tell anything about how much the corresponding 
gauge configurations are differing from each other 
(irrespective of a global relative gauge transfor- 
mation) we would like to interprete our finding 
of shrinking distributions as a weakening of the 
Gribov problem with increasing 'physical' lattice 
size. 

Moreover, we have plotted the relative devia- 
tion 



SD L ( Pmm ) = (D^l - D&Mok (11) 

for the gluon propagator (see Fig. and anal- 
ogously for the ghost propagator (see l.h.s. of 
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FIG. 7: Distributions of the deviation of the gauge 
functional values for different gauge copies relative 
to the best copy per configuration. FOR method at 
P = 2.40 for lattice sizes 8 4 and 16 4 . 
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FIG. 8: Distributions of ghost propagator values 
at lowest non-trivial momentum for different gauge 
copies as in Fig. |7] 



Fig. I1UII as a function of the inverse linear lat- 
tice size 1/L , both determined at the minimal 
momentum p m in- Here we have used data for 
fixed (3 = 2.4 and lattice sizes from L = 5 up 
to L = 20. In close correspondence to our ob- 
servations presented in Figs. |21and|5]we see that 
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p=2.4 



1/L 



FIG. 9: 



D 



(be) 



Relative deviation 5D_l(p„ 

-.(6c) 



(fc) _ 
SOR 

] i)r)/^for m percent for the gluon propagator 
D for various linear lattice sizes L and smallest 



non-vanishing momentum p„ 
(/3 = 2.4). 



= (2/a) sin(7r/L) 



the Gribov copy effect becomes weaker (stronger) 
for increasing (decreasing) 'physical' lattice size 
and correspondingly decreasing (increasing) min- 
imal momentum, at least up to a certain value of 
the lattice size (^15). One would of course need 
larger values of L to make a reliable conclusion 
about the limit L — * oo. Anyway, at our largest 
lattice value L = 20 the Gribov copy effect is still 
quite strong. 

For the ghost propagator, where the signal to 
noise ratio is more favourable, we have found an 
analogous behavior also for the multiple on-axis 
momenta k = 2, 3, 4 (see r.h.s. of Fig. I1L)[I . 

In [13| two of us have reported on rare Monte 
Carlo events with exceptionally large values of 
the ghost propagator occuring for the SOR gauge 
fixing method for larger values. In Fig. 1111 we 
show some time histories for the gluon and ghost 
propagators for (3 = 2.5 and a 16 4 lattice, com- 
paring be SOR with be FOR. We see that for 
the 'best copy - flips on' case (FOR) the fluctu- 
ations for both propagators are smaller. But for 
the ghost propagator the effect of exceptionally 
large values, in general related to small eigenval- 
ues of the F-P operator is still there. 

Concluding we show the form factors of the 
gluon propagator p 2 D(p) and of the ghost prop- 
agator p 2 G(p) in physical units as a function of 
the physical momentum for fixed /3 = 2.4 and 
lattice sizes varying from 10 4 to 20 4 . We have 



rescaled the gluon propagator values D(p) with 
factors a 2 and g$ and the ghost propagator 
G(p) with a 2 , respectively, in order to translate 
to the corresponding continuum (bare) propaga- 
tors (compare with [22]). To estimate the lattice 
spacing in physical units we have used the string 
tension: a 2 a = .071 |2^ with the standard value 
yj~d = 440MeV. The form factor results for both 
methods be SOR and be FOR are shown to- 
gether in Fig. [U A gain the figure shows clear 
Gribov copy effects for both the propagators and 
not only for the ghost propagator. We did not ap- 
ply any overall renormalization here. The statis- 
tics collected for these runs is listed in Table Hill 



FOR 


L 


5 


6 


8 


9 


10 


12 


15 


16 


18 


20 


N cori f 


1000 


1000 


800 


600 


600 


500 


400 


356 


200 


200 


SOR 


L 


5 


6 


8 


9 


10 


12 


15 


16 


18 


20 


Nconf 


1000 


1000 


800 


500 


500 


400 


400 


370 


100 


100 



TABLE III: Statistics for the measurements at dif- 
ferent L and f3 = 2.4. 



V. CONCLUSIONS 

In this paper we have demonstrated that there 
is a visible Gribov problem for the ghost prop- 
agator as well as for the gluon propagator com- 
puted in SU(2) lattice gauge theory within the 
Landau gauge. In order to show this we have en- 
larged the gauge orbits of given Monte Carlo gen- 
erated gauge fields by non-periodic Z(2) transfor- 
mations, flipping all links in a given direction on 
a slice orthogonal to that. This allows a precon- 
ditioning which maximizes the gauge functional 
before applying the overrelaxation algorithm. 

We have found indications for a weakening of 
the Gribov copy effect both going to larger mo- 
menta at fixed volume and also increasing the 
lattice size L while correspondingly lowering the 
minimal non-zero momentum, at least up to a 
certain value of the lattice size (^$15). However, 
one would need larger values of L to draw a reli- 
able conclusion about the limit L — > oo. 

We have not shown the momentum scheme 
running coupling which can be determined from 
the form factors of the propagators discussed 
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FIG. 10: Relative deviation 8Gl(p) = (Gg^p — GpQ R )/GpQ R in percent for the ghost propagator G at 
/3 = 2.4 for various linear lattice sizes L and the smallest non-vanishing momentum p m in = (2/a) sin(7r/L) 
(left) as well as for on-axis momenta p(k), k — 2,3,4 (right). 
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FIG. 11: Time histories for D {bc) {p min ) (left) and G {bc) (p min ) (right) for both SOR and FOR methods at 
(3 = 2.5 and 16 4 lattice. 



here assuming that the renormalization factor for 
the ghost-gluon vertex is constant. This will be 
discussed in a future paper, where we want to 
present data for larger lattices and a larger spec- 
trum of (off-axis) momenta. 
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